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Abstract 

We consider the most general four fermion operators in QCD for two and three 
massless flavors and study their renormalization in the vicinity of the Fermi 
surface. We show that, asymptotically, the largest coupling corresponds to 
scalar diquark condensation. Asymptotically the direct and iterated (molec- 
ular) instanton interactions become equal. We provide simple arguments for 
the form of the operators that diagonalize the evolution equations. Some solu- 
tions of the flow equations exhibit instabilities arising out of purely repulsive 
interactions. 
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1. Recently there has been renewed interest in the problem of hadronic matter at high 
baryon density. It was realized early on that at very high density, asymptotic freedom and 
the presence of a Fermi surface implies that attractive interactions between quarks cause 
a BCS instability, and cold quark matter is a color superconductor [I]. More recently, the 
problem was studied again, and it was emphasized that non-perturbative effects, instantons, 
could lead to gaps on the order of 100 MeV PJ| . A particularly interesting case is QCD with 
three light flavors. In this case the preferred order parameter involves a coupling between 
color and flavor degrees of freedom [|J . The resulting coherence leads to large gaps even for 
weak interactions. In addition to that, color- flavor locking implies that chiral symmetry is 
broken. The usual quark condensate is non-vanishing even in the high density phase, as is 
a gauge- invariant AB = 2 condensate, signalling true superfluidity ||||. 

Most of these calculations were performed in the mean field approximation. (In the 
instanton model, a number of refinements have been included, see [§,0.) One assumes a 
specific form of the interaction (instantons, one gluon exchange, . . .), makes an ansatz for 
the form of the condensate, and solves the gap equation in some approximation. Clearly, it is 
desirable to analyze the structure of the theory for the most general form of the interaction, 
and to have a guiding principle for constructing a systematic expansion scheme. 

In this context, the renormalization group approach to cold Fermi systems appears par- 
ticularly promising f^,§ . In a cold Fermi system, the only relevant interactions take place in 
vicinity of the Fermi surface. The corresponding excitations are quasi-particles and quasi- 
holes, described by an effective action of the form 

S eff = Jdtd 3 p^(idt-«p)-e F ))^. (1) 

We can analyze the general structure of the quasi-particle interactions by studying the 
evolution of the corresponding operators as we successfully integrate out modes closer and 
closer to the Fermi surface. The main result of this analysis is that, in general, four fermion, 
six fermion, and higher order interactions are suppressed as we approach the Fermi surface. 
This fixed point corresponds to Landau liquid theory ||. The only exception is the four- 
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fermion operator that corresponds to two particles from opposite corners of the Fermi surface 
(Pf, ~Pf) scattering into (p' F , —p' F ). This kind of scattering leads to a logarithmic growth of 
the coupling constant as we approach the Fermi surface, and the well known BCS instability. 

2. For QCD, this approach was first employed by Evans, Hsu, and Schwetz ||10|| . These 
authors concentrated on the one gluon exchange interaction, but as we will see in a moment, 
their results are sufficient to deal with the most general quark-quark interaction for QCD 
with three (or more) flavors. We should emphasize here that we will restrict ourselves to 
massless QCD, a spherical Fermi surface, and local operators invariant under the appropriate 
chiral symmetry. 

Following |1(J, we take the basic four fermion operators in Nf = 3 QCD to be 

oI l = ^li^l)\ o° lr = tyLiMtyiaM (2) 

Each of these operators comes in two color structures, for example color symmetric and color 
anti-symmetric 

(r^ b )(^ d )(6 ab 5 cd ±5 ad 5 bc ). (3) 

Nothing essentially new emerges from considering superficially different isospin structures, 
or different Dirac matrices. All such structures can be reduced to linear combinations of 
the basic ones (§), or their parity conjugates, by Fierz rearrangements. In total, we have to 
consider eight operators. 

These operators are renormalized by quark-quark scattering in the vicinity of the Fermi 
surface. This means that both incoming and outgoing quarks have momenta pi,f>2 — ±p 
and P3,Pa — ±<f with \p\, |g| ~ pp. We can take the external frequency to be zero. A graph 



with vertices Ti and T2 then gives [1C 
G\G 2 I [Ti)i'i{^i)k'k 



(7o)v(To)« - ^(l)ij(l)ki 



(r 2 ),y(r 2 V (4) 



with I = -^ji 2 log(A/#/ Kjjv)- Here [Am, Ajjy] is the range of momenta that was integrated 
out. We will denote the density of states on the Fermi surface by iV = /i 2 /(27r 2 ) and the 



logarithm of the scale t = log(Ajji/ Ajjv) as in JT0|. The renormalization group does not mix 



LL and LR operators, and it also does not mix different color structures. This means that 
the evolution equations contain at most 2x2 blocks. For completeness, we reproduce the 
results of [fj 

d{G$ L + G^ L ) N , nLL n LL\2 tr\ 

- dt ~ = + G * } (5) 

d(G LL - 3GP) = _ n{g ll _ 3G LLy (6) 



dt 

d(G% R + 3G^ R 



(7) 



dt 

r<LR\ o AT 

-(G^-O 2 (8) 



d(G% R - Gf R ) 2N 



dt 

In this basis the evolution equations are already diagonal. The coupling G\ = G LL + G\ L 
evolves as 

Gl{t) = 1 + (N/^GMt (9) 
with analogous results for the other operators. Note that the evolution starts at t — and 
moves towards the Fermi surface t — oo. If the coupling is attractive at the matching scale, 
Gi(0) > 0, it will grow during the evolution, and reach a Landau pole at t c = 3/(NGi(0)). 
The location of the pole is controlled by the initial value of the coupling and the coefficient in 
the evolution equation. If the initial coupling is negative, the coupling decreases during the 
evolution. The second operator in (f|) has the largest coefficient and will reach the Landau 
pole first, unless the initial value is very small or negative. In this case, subdominant 
operators may dominate the pairing. 

The particular form of the operators that diagonalize the evolution equations is easily 
understood. Let us first order the operators according to the size of the coefficient in the 
evolution equations 

O dom = {^lIo^l) 2 - {^lI^l) 2 , (10) 
O sub)1 = {^l1o^l){^r1o^r) - ^(4>lj4>l)(4>ri4>r), (11) 
O sub , 2 = {^lIo^Pl) 2 + \{^l1^l) 2 , (12) 



O mar = {^LlO^l) {^RIo^r) + {^LI^l) {^RI^r) ■ (13) 

This result can be made more transparent by Fierz rearranging the operators. We find 



Odom 




(14) 


O su b,l 




(15) 


O su b,2 


= ^ty L CEil> L )$ L CE$ L ), 


(16) 


O m ar 


= ^{i'LC'joipR){i!RC'yoi! L ) + .... 


(17) 



This demonstrates that the linear combinations in (|T0l - |l3"|) correspond to simple structures 
in the quark-quark channel. It also means that it might have been more natural to perform 
the whole calculation directly in a basis of diquark operators. 

The full structure of the (LR) operators is S ub,i,O ma r = (i ) C r yrs,A' l l>)('^'CjTs,A'^) + 
(iiC , y'y5TA,s' i P)(4>Cjj5TA,s4>) where ts,a are symmetric/anti-symmetric isospin generators. 
Note that because the two structures have different flavor symmetry, the flavor structure 
cannot be factored out. 

The dominant operator corresponds to pairing in the scalar diquark channel, while the 
subdominant operators contain vector diquarks. Note that from the evolution equation 
alone we cannot decide what the preferred color channel is. To decide this question, we 
must invoke the fact that "reasonable" interactions, like one gluon exchange, are attractive 
in the color anti-symmetric repulsive in the color symmetric channel. Indeed, it is the color 
anti-symmetric configuration that minimizes the total color flux emanating from the quark 
pair. If the color wave function is anti-symmetric, the dominant operator fixes the isospin 
wave function to be anti-symmetric as well. 

The dominant operator does not distinguish between scalar and pseudoscalar diquarks. 
Indeed, for Nf > 3 the basic 4-quark operators consistent with chiral symmetry also ex- 
hibit an accidental axial baryon symmetry, under which scalar and pseduscalar diquarks are 
equivalent. For Nf = 3 this degeneracy is lifted by (formally irrelevant) six-fermion opera- 
tors |0] . The form of the dominant operator indicates the existence of potential instabilities, 
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but does not itself indicate how they are resolved. For example, to see whether color-flavor 
locking Q is important we should use the operator as input to a variational calculation. 

3. For two flavors, we have to take into account additional operators. At first hearing it 
might seem odd that with fewer basic entities we encounter more basic operators. It occurs 
because for Nf = 2, but not for larger values, two quarks of the same chirality can form a 
chiral SU(2) x SU(2) singlet. Related to this, for Nf = 2 we have additional U(1)a violating 
four fermion operators. These operators are induced by instantons. 

For Nf = 3, instantons are six fermion operators, and they are irrelevant in the technical 
sense. This does not mean they are physically irrelevant, particularly since they break 
a residual symmetry. The formation of the gap will cause the evolution of the couplings 
to stop, and the instanton coupling remains at a finite value. Instantons have important 
physical effects, even for Nf = 3 0. Most notably, instantons cause quark-antiquark pairs 
to condense (even in the high density phase), and lift the degeneracy between the scalar and 
pseudoscalar diquark condensates. 

The new operators are 

O s = det(^L), T = det(^S^) (18) 

Both operators are determinants in flavor space. For quark-quark scattering, this implies 
that the two quarks have to have different flavors. The fact that the flavor structure is fixed 
implies that the color structure is fixed, too. For a given (qq) spin, only one of the two 
color structures contributes. Finally, both quarks have to have the same chirality, and the 
chirality is flipped by the interaction. 

These considerations determine the structure of the evolution equations (see Fig. 1). Two 
left handed quarks can interact via one of the instanton operators, become right handed, 
and then rescatter through an anti-instanton, or through one of the U(1)a symmetric RR 
operators. The result will be a renormalization of the LL vertex in the first case, and 
a renormalization of the instanton in the second. The iterated instanton-anti-instanton 



interaction was discussed in great detail in [pjj| , and was argued to play an important role 



in the high temperature |TTJ and high density phase of QCD 0. We should note that the 
flavor structure will always remain a determinant. Even though instantons generate all the 
Dirac structures in (0), the color-flavor structure is more restricted. 

Evidently, instantons do not affect the evolution of the LR couplings at all. The evolution 
equations of the LL couplings are modified to become (henceforth we drop the subscript LL 
from) : 

dG Q N ' 



dt 2 



- G 2 + 2G Gi - 5G 2 -K 2 S + 2K S K T - 5K 2 (19) 



2 (-Go + Gi) K s + 2 (Go - 5G,) K T (21) 



d -^ = *\ l -Gl- ^GoG, + l M + -Kl - ™K S K T + H**) (20) 
dt 2\3 3 3 i 3 5 3 3 T J V ' 

dK s _ N 
~df ~ ~2 v 

= y || (G - 5G.) K 5 + | (-5G + 13G,) K T j (22) 

These equations can be decoupled as 

rfgi _ JV 

dK x 2N 



dt 3 
dG 2 

dt 
dK 2 



G\ + X?) , (23) 
GxKx, (24) 



-iV (Gl + Kl) , (25) 
-2N G 2 K 2 , (26) 



where G\ = Go + Gj, K\ = Ks + Kt and G 2 = Go — 3Gj, = ^Cs — 3Kt- The equations 
for G, K decouple even further. We have 

d{G 2 + K 2 ) 



dt 

d(G 2 - K 2 ) 



-N (G 2 + K 2 y (27) 
-N(G 2 -K 2 ) 2 , (28) 



dt 

as well as the analogous equation for Gi,K\. These differential equations are now trivial to 
solve, leading to 

G ^-l{^TWt + dm}^ W 

K M = Ut^-tJ™)> (3°) 



2 {a + Nt b + Nt 
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again with the analogous result holding for Gi,K 1 . Here, a,b = (G 2 (0) ± K 2 (0))~ 1 . The 
result implies that G 2 and K 2 will grow and eventually reach a Landau pole if either a or b 
is positive. The location of the pole is determined by the smaller of the values, t c = —a/N 
or t c = —b/N. The same is true for G\ and K\, but the couplings evolve more slowly, and 
the Landau pole is reached later. 

At this level a number of qualitatively different scenarios are possible, depending on the 
sign and relative magnitude of G(0) and K(0), see Fig. 2 (henceforth we drop all subscripts). 
If G(0) and K(0) are both positive then they will both grow, and the location of the nearest 
Landau pole is determined by G(0) +K(0). The asymptotic ratio of the two couplings is 1. 
If G(Q) and K(Q) are both negative, and the magnitude of G(0) is bigger than the magnitude 
of K(0), then the evolution drives both couplings to zero. These are the standard cases. 
Attraction leads to an instability, and repulsive forces are suppressed. 

More interesting cases arise when the sign of the two couplings is different. The case 
G(0),K(0) < and \K(0)\ > |G(0)| is especially weird. Both G{0),K{0) are repulsive, but 
the evolution drives G(0) to positive values. Both couplings reach a Landau pole, and near 
the pole their asymptotic ratio approaches minus one. Similarly, we can have a negative 
G(0) and positive K(0) with K(0) > |G(0)|. Again, the evolution will drive G(0) to positive 
values. 

4. The dominant and sub-dominant instanton operators are 



Odorn = det 
S ub = det 



Upon Fierz rearrangement, we find 



Ml) 2 - (4>r^l) 2 \ , (31) 
(Ml) 2 + 1$i£iI>l) 2 ]. (32) 



O dom = 2(i/> l Ct 2 iI> l )$ r Ct 2 $ r ), (33) 
O sub = ^ l Ct 2 ^ l )(4, r Ct 2 ^ r ), (34) 

corresponding to scalar and tensor diquarks. Both operators are flavor singlet. Overall 
symmetry then fixes the color wave functions, anti-symmetric 3 for the scalar, and symmetric 



6 for the tensor. The dominant pairing induced by instantons is in the scalar diquark channel, 
the only other attractive channel is the tensor. All this neatly confirms the scenario discussed 
in Note that condensation in the tensor channel violates rotational symmetry. As a 

result, the gap equation has additional suppression factors and the gap is very small |2|. 

Just as we found for Nf > 3, there is an appealing heuristic understanding for the 
amazingly simple behavior of the evolution equations, obtained by focussing on the diquark 
channels. Instantons distinguish between scalar diquarks with positive and negative parity. 
G + K corresponds to the positive parity operator (ipC^ip) and G — K to the negative parity 
(i/jCi/j). The asymptotic approach of G/K — > 1, then corresponds to the fact that scalar 
diquark condensation is favored over pseudoscalar diquark condensation. This is always 
the case if K(0) > 0. We also understand the strange case G(0),K(0) < and |i^(0)| > 
|G(0)|. In this case the interaction for scalar diquarks is repulsive, but the interaction in 
the pseudoscalar channel is attractive and leads to an instability. Note that this can only 
happen if we have the "wrong" sign of the instanton interaction, i.e. for 9 = n. Similarly, we 
can understand why the asymptotic ratio of the molecular (instanton-anti-instanton) and 
direct instanton couplings approaches G/K = ±1. Instantons induce a repulsive interaction 
for pseudoscalar diquarks. During the evolution, this coupling will be suppressed, whereas 
the attractive scalar interaction grows. But this means that in the pseudoscalar channel, 
the repulsive (instanton) and attractive (molecular) forces have to cancel in the asymptotic 
limit, so the effective couplings become equal. 

We have not made an attempt to match all the coupling constants to a realistic model at 
the UV scale. To do so would require an understanding of the instanton density, the relevant 
value of a s , the screening mechanism, and many other things. From the form of the instanton 



vertex we can fix the ratio of the two instanton-like couplings, Kt/K$ = 1/(2N C — 1) flT2| , |T3" 
This again shows that the tensor channel is not expected to be important. Also, both one 
gluon exchange and higher order instanton effects give Gq L > 0, —G\ L > 0. Thus the 
favored scenario is that both instanton and U(1)a symmetric couplings flow at the same 
rate, and pairing is dominated by scalar diquarks. 



5. In summary, we find that the renormalization group analysis broadly supports the 
findings of 0,0]. The dominant coupling corresponds to scalar diquark condensation, the 
sub-dominant coupling to tensor diquarks. Asymptotically, U(1)a breaking and U(1)a sym- 
metric couplings flow at the same rate. Since instantons have a definite isospin structure, 
there is one flavor symmetric coupling that evolves independently of instantons. Asymp- 
totically, this coupling also flows at the same rate. In principle there is the possibility of 
flavor and color symmetric diquark condensates, but none of the model interactions so far 
proposed is attractive in that channel. 

The analysis presented here is incomplete in several ways. Realistic interactions (one- 
gluon exchange, instantons, . . .) are momentum dependent, and that should be included in 
the evolution. This complication is particularly important for one-gluon exchange, because 
in perturbation theory the diagram is not completely screened, and has a divergence for 
small momentum transfers. For Nf = 3 a self-consistent calculation ought to be possible, 
because color-flavor locking completely screens the interaction. 

In any case, the renormalization group only determines the running of the couplings from 
the matching point down to some infrared scale. If a subdominant coupling is unusually 
large at the matching scale, it might still dominate the pairing. To decide the form of the 
pairing and determine the gap, given the couplings, still requires a variational calculation 
along the lines of [|],f|. 

Acknowledgments: After this work was begun, we learned that the authors of |10[ have 
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by NSF-PHY-95 13835. 
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FIGURES 



FIG. 1. Chiral structure of the evolution equations. 

FIG. 2. Solutions to the evolution equations for cases with different relative strength of the 
U(1)a violating interaction K (dashed line) and U(1)a conserving interaction G (solid line). 
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